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^ ! Abstract 

Q . 

We present a first physical application of Tsallis' generalized entropy to the ther- 

pg ' modynamics of self-gravitating systems. The stellar system confined in a spherical 

^ ■ cavity of radius rg exhibits an instability, so-called gravothermal catastrophe, which 

^J , has been originally investigated by Antonov (Vest. Leningrad Gros.Univ. 7 (1962) 

^ ; 135) and Lynden-Beh & Wood (Mon.Not.R.Astron.Soc. 138 (1968) 495) on the 

t^ ■ basis of the maximum entropy principle for the phase-space distribution function. 

In contrast to previous analyses using the Boltzmann-Gibbs entropy, we apply the 
Tsallis-type generalized entropy to seek the equilibrium criteria. Then the distri- 
bution function of Vlassov-Poisson system can be reduced to the stellar polytrope 
system. Evaluating the second variation of Tsallis entropy and solving the zero 
eigenvalue problem explicitly, we find that the gravothermal instability appears in 
'^ ' cases with polytrope index n > 5. The critical point characterizing the onset of 

^ ■ instability are obtained, which exactly matches with the results derived from the 

Q I standard turning-point analysis. The results give an important suggestion that the 

Tsallis' generalized entropy is indeed applicable and viable to the long-range nature 
of the self-gravitating system. 
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1 Introduction 



In any subject of astrophysics and cosmology, many-body gravitating systems 
play an essential role. In general, dynamics of such systems are quite difficult to 
understand and the long-range nature of gravitational interaction prevents us 
from applying the statistical mechanics. Thermodynamics of self-gravitating 
systems also shows some peculiar features such as a negative specific heat 
and an absence of global entropy maxima, which greatly differ from usual 
thermodynamic systems. 

To see the peculiarity of self-gravitating systems, consider a system confined 
within a spherical adiabatic wall. We assume that the particles in this sys- 
tem interact via Newton gravity and bounce elastically from the wall. Under 
keeping the energy and the total mass of the system constant, a thermody- 
namic description of self-gravitating systems leads to an interesting conclusion. 
When the central mass density is sufficiently high, no equilibrium state ex- 
ists and the system can persistently undertake a strong central condensation. 
This instability is known as gravothermal catastrophe, originally investigated 
by Antonov[l] and Lynden-Bell and Wood[2] (see also Refs.[3-6]). They define 
the entropy of distribution function in phase-space and seek the equilibrium 
condition whether the system exhibits the maximum entropy state. In their 
analyses, they treat the Boltzmann-Gibbs entropy. 



5bg = - / / In / d^xd^p, 



where f{x,p) is the distribution function in phase-space. 

While the gravothermal catastrophe found in early studies has been widely 
accepted as a fundamental astrophysical process and plays an essential role in 
the dynamics of globular clusters [7-10], a naive but natural question arises. 
Why should we use the Boltzmann-Gibbs entropy when looking for a prob- 
able entropy state ? The choice of entropy severely restricts the functional 
form of the distribution function /. In a spherically symmetric system with 
isotropic velocity distribution, the entropy (1) leads to the isothermal gas dis- 
tribution. The equilibrium configuration, however, cannot be unique in the 
self-gravitating systems. There remains a possibility of another choice of en- 
tropy in order to determine the most probable state. 

Indeed, in the view of statistical mechanics and thermodynamics, it has re- 
cently been known that the standard formalism based on the Boltzmann- 
Gibbs entropy cannot deal with a variety of interesting physical problems and 
serious difficulties arise when applying the Boltzmann-Gibbs statistical me- 
chanics (e.g, [11-13]). After introducing a family of generalized entropies by 
Tsallis [14], a new framework of thermodynamic structure has been exten- 



sively discussed [15-19] and some physical applications has been presented 
successfully (e.g, [20-22]; see [23,24] for comprehensive reviews). Although 
there still remains a fundamental issue on the consistency with the thermo- 
dynamic relations [19], the new formalism is expected to be applied in a quite 
wide area including physics, astronomy, biology, economics, etc. Especially, in 
many astrophysical problems involving the long-range nature of gravity, the 
nonextensive generalization of statistical and/or thermodynamic treatment 
should deserve further consideration. At present, however, no specific exam- 
ples have been reported theoretically, except for the experimental result [25] 
or observational evidence [26]. 

If one applies Tsallis' generalized entropy to the above problem in self-gravitating 
systems, the equilibrium condition is explored by using 



^Jif- /) d'^d'P^ (2) 



s. 



instead of (1). Here, the parameter q is chosen as g 7^ 1, as a possible general- 
ization of the Boltzmann-Gibbs entropy. Using (2), an attempt to determine 
the most probable state has been made by Plastino and Plastino [27,28]. They 
found that the equilibrium configuration reduces to the polytropic gaseous sys- 
tem, which has been widely utilized in a study of the stellar structure [29,30]. 

Based on their result, we then pursue to investigate the stability of self- 
gravitating system confined within a box and try to answer the crucial ques- 
tion; how the equilibrium condition is altered when applying the entropy (2) ? 
More basically, does the Tsallis type entropy correctly predict the equilibrium 
condition of a stellar polytrope system ? 

This paper is organized as follows. In section 2, we briefiy review that the 
most probable state determined by the Tsallis entropy reduces to the stellar 
polytrope. Then we proceed to the stability analysis based on the standard 
turning-point analysis in section 3. In section 4, stability/instability criterion 
is re-derived from the second variation of entropy. The stability/instability 
criterion can be obtained by solving zero-eigenvalue problem, which exactly 
matches the criterion from the standard turning-point analysis. Also, the geo- 
metrical construction of stability/instability criterion is discussed in terms of 
homology invariants. Finally, section 5 is devoted to conclusions and discus- 
sion. 



2 Tsallis entropy and stellar polytrope 



In this section, for the notational convenience and the subsequent analyses, 
we start to review the work by Plastino & Plastino [27] that the extremum 
state of the Tsallis entropy just reduces to the stellar polytrope in section 2.1. 
Adopting this extremum state, the equilibrium configuration is then deter- 
mined by solving the conditions of hydrostatic equilibrium in section 2.2. A 
family of equilibrium sequence referred to as the Emden solutions is obtained 
and characterized by the homology invariants, which is subsequently used in 
the stability analysis. 



2. 1 The principle of maximum entropy 



Consider a system containing A^ particles which are confined within a cavity of 
hard sphere. The radius of cavity is r^ and each particle is assumed to have the 
same mass m. Then, all the information of this system can be described by 
a A^-body distribution function /at (aji, ■ ■ ■ , Xp^; Pi, ■ ■ ■ , Pj^; t), defined in the 
6A^-dimensional phase-space. The evolution of such a system is governed by 
the collisionless Boltzmann equation, which is generally intractable. Instead of 
using the full distribution, it is better to employ a coarse-grained distribution 
function whose value at any phase-space point {x,p) is averaged over some 
specified volume centered on {x,p). In this treatment, the correlation between 
particles is smeared out. Thus, the system can be regarded as a gravitating 
gaseous system. The phase space distribution is simply described by the "one- 
particle distribution function" f{x,p;t), which greatly reduces the problem 
to a tractable level. 

Adopting the mean-field treatment, we now investigate the equilibrium config- 
uration and stability of a gravitating gaseous system with a help of statistical 
mechanics. In this approach, the equilibrium configuration can be determined 
by the principle of maximum entropy, which specifies the distribution function 
/ that maximizes the entropy. As is well-known, however, the global maximum 
does not exist in self-gravitating systems [1,2] [8]. Hence, we try to look for the 
local extrema 6S = and seek the criterion whether these extrema are really 
local maxima or not. Previously, most previous work has extensively discussed 
this issue on the basis of the Boltzmann-Gibbs entropy (1), in which case the 
extremum solution reduces to the isothermal gas system. In this paper, with 
a particular attention to the Tsallis- type generalized entropy (2), we seek the 
equilibrium configuration under the mass and the energy conservation: 



M 



I f d^xd^v, (3) 



E = K + U = ^ fv^f dJ^xd?v + ]-[ $(x) / d^xd^v, (4) 



where the quantity $ denotes gravitational potential: 



$(a;) = -G / p^' rf^yrf^-y, (5) 



with V = p/m. Here, we adopted the standard definition of mean value for the 
mass M and the energy E. The remark on the use of other definitions such as 
g-generalized mean value [17] is discussed in section 5. 

The extremum entropy state can be derived by varying Sq with respect to /. 
Using the Lagrange multipliers a and /9, the extremum solution subject to the 
constraints (3) and (4) is obtained from 

5Sq-a5M -P5E = Q, (6) 

which leads to 



/ {-^(^^'~' - 1) - a - /3 i^-v^ + $) I 5/ (fxd'v = 0. 



(7) 



Here the relation J S^ f d^xd^v = J ^6f d^xd^v is used in deriving the above 
expression. Since the constraint (6) must be satisfied independently of the 
choice of 6f, we obtain 

1 .(^qP^^^l) + a + (3(lv' + <^)=0, (8) 



g-r^" ' ' \2 



which reduces to the following distribution function: 



f{x,v) = A 



% - ^{x) - -v^ 



(9) 



where we define the constants A and $o- 



The functional form of the distribution function (9) implies that the extremum 
solution is indeed equivalent to the polytropic gaseous systems [8] [27]. The 
density profile p(r) at radius r = |a;| can be expressed using (9): 



p{x) = J fd^v, 

= 4v^7r B (^, -^) A {$0 - $(a;)}i^+i (11) 

with B{a, b) being the P function. On the other hand, in the case of isotropic 
velocity distribution, the pressure becomes 

P(x) ^ j\v' f d'v, 

= (^ + ^) /^(^) i^^o-Hx)}. (12) 

Thus, these two equations lead to the relation 

P OC p(55-3)/(39-l) , (13) 

which corresponds to the polytropic equation of state, P oc p^+^/", and the 
polytrope index n is connected with Tsallis' g-parameter as follows: 

n=^ + |. (14) 

q — 1 2 

Note that the polytrope index n should be positive in any astrophysical system. 
It has also been argued that the values of n in the interval < n < 3/2 
are unphysical[8], which restricts the parameter q to be larger than unity. 
Further, a simple argument shows that the locally maximum entropy can only 
be attained in the spherically symmetric configuration [1,2]. Although the 
extremum solution (9) does not restrict any symmetry and any value of the 
parameter q, we hereafter restrict our attention to the spherically symmetric 
case with g > 1. 

The resultant distribution (9) and the relation (14) further imply that the 
polytrope gas is equivalent to the isothermal gaseous sphere in the limit g — > 1 
or n — > +00, which is in fact obtained in adopting the Boltzmann-Gibbs 
entropy (1). To check this consistency, let us take the limit g ^ 1. Introducing 
the new constant T: 

T^iq-l)\AV2nB -,-L-] a\ , (15) 



one can rewrite the distribution function (9) using the relation (14) and (15) 
as 



f(r, v) = — = 

^ ^ 4^2715(3/2, n- 1/2) 



X 



p(r) [^ vy2 T"-'/' 



{(n-3/2)Tpi/n(r)}'/' 



1 - 



(n- 3/2) Tpi/«(r) 



• (16) 



Also the polytropic equation of state becomes 



Pir) = (^^^) Tp^^^'-{r). (17) 



Then using the fact that 5(3/2, n — 1/2) — > {ix/Arv'Y^'^ in the hmit n -^ oo, 
the distribution function asymptotically approaches 



/(r,.)"=±?=p^p(r)e-'/-, (18) 

and the equation of state reduces to that of the isothermal gas, P = pT. The 
expression (18) shows that the velocity distribution is indeed Maxwelhan and 
functional form of the distribution is fully specified by the 'temperature', T 
(velocity dispersion). 



2.2 Emden solution 



While the extremum state is determined by the variation of entropy, the dis- 
tribution function (9) does not yet completely specify the equilibrium config- 
uration. In the expression (9), the gravitational potential $ appears, which 
implicitly involves / again (see eq.(5)). To break this roundabout, an explicit 
expression for $ or p is needed. 

In the spherically symmetric configuration, the gravitational potential given 
by (5) satisfies the following Poisson equation: 

r'^ dr \ ar I 



Combining (19) with (11), we obtain the ordinary differential equation for $. 
Alternatively, a set of equations are derived using (11), (12) and (19): 



^^-^,W. (20, 

^ = 4.,M.^ (21) 



which represent the hydrostatic equihbrium. The quantity m{r) denotes the 
mass evaluated at the radius r inside the wall. We then introduce the dimen- 
sionless quantities: 



which yields the following ordinary differential equation: 

e" + jO' + ^" = 0, (23) 



where prime denotes the derivative with respect to ^. The quantities pc and 
Pc in (22) are the density and the pressure at r = 0, respectively. To obtain 
the physically relevant solutions, we solve equation (23) with the following 
boundary condition: 

^(0) = 1, ^'(0) = 0. (24) 



A family of solutions satisfying (24) is referred to as the Emden solution^ which 
is well-known in the subject of stellar structure and details of the solutions 
can be found in standard textbooks (e.g., see Chap. IV of ref. [29]). Except for 
the few cases n = 0, 1 and 5, the Emden solution cannot be expressed in terms 
of the elementary functions. So the solution with general index n is obtained 
numerically. 

Notice the fact that equation (23) is invariant under the transformation ^ — > 
A^^ 9 -^ A^'^^^^^^^O, where A is an arbitrary constant. This implies that, given 
a solution with some value of ^(0), we can obtain the solution with any other 
value of 6{0) by simple rescaling. Therefore, only one of the two integration 
constants in (23) is actually non-trivial. This fact allows us to reduce the 
degree of equation from two to one by suitable choice of variables. One such 
set of variables is 



_ 4nr'p{r) _ ^9" 
m[r) u 

F[r) r u 

which are called homology invariants[29, 30]. That is, for a fixed n, all the 
solutions depicted in [u, i))-plane lie on the same trajectory. In terms of these 



variables, the equation (23) can be written as 

u dv (n + l)(-u — 1) + f 
V du (n + 1)(3 — u) — nv 



(27) 



Figure 1 shows the Emden solutions for various polytrope indices. As indicated 
by the boundary condition, (24), all the trajectory start from (m, f) = (3,0), 
which represents the center of the configuration, r = 0. As the radius r 
increases, the trajectory monotonically moves to the upper-left direction in 
(m, w)-plane, as marked by arrow. For larger radius, the solution with poly- 
trope index n > 5 spirals around a fixed point, while the trajectories with 
n < 5 continue to approach the point {u,v) = (0, oo). The marginal case 
is n = 5. The trajectory monotonically changes and asymptotically reaches 
{u,v) = (0,6). These differences can be explained by the behavior of outer 
envelope p(r) as follows. In cases of n < 5, the density p{r) falls off rapidly 
and it eventually vanishes at a finite size, where the pressure also becomes 
zero. Since the mass is finite and the polytropic relation P oc p^+^/" holds, 
the quantities {u,v) becomes {u,v) = (0, oo) from (25) and (26). On the other 
hand, when n > 5, the outer profile of the density falls off less steeply than 
p oc r~^[8]. In this case, the trajectory in (■u,f)-plane cannot reach u = 0. 
It must be bounded within the interval < -u < 3 and shows the oscillatory 
behavior, as depicted in Figure 1. 

The characteristic feature seen in Figure 1 provides us an important suggestion 
about stability of the polytrope gas sphere. This point will be clarified by the 
turning-point analysis and the second variation of entropy in the subsequent 
section. 



3 Stability/instability criterion from the turning-point analysis 



We are specifically concerned with the stability/instability of static equilibria 
from the thermodynamic point of view. For this purpose, following the dis- 
cussion in [2,5,6], we first apply the standard-turning point analysis to the 
equilibrium sequence obtained in the previous section. 

Recall that the system is confined within the spherical adiabatic wall of radius 
Te. For fixing polytrope index n, the local extremum satisfying 6Sq = is 
characterized by the energy E, mass Mand radius of rigid sphere rg. In the 
absence of instability, any value of E[—oo < E < oo), M(0 < M < oo) and 
re(0 < Te < oo) are accommodated by a suitable choice of pc and Pc, but, 
there exists a lower bound on the dimensionless quantity Eve/GM'^ for the 
polytrope gas solutions. To show this, we compute the total energy contained 



within the hard sphere of a radius rg. The kinetic and potential energies, K 
and U are respectively expressed as 



K = - drAnr'^ P{r) 





(2J 



and 



U-- 



Gm dm 


G f dr d 
2 J r dr 




r dr 



dr 



A-nriPe - 3 / drAnr'Pir). 



[m 



GM^ 



dr 



Gm^ 



(29) 



Combining these results and using (57) in appendix A, the total energy E can 
be written in a compact form: 



E = K + U = A-wrlPe -^ I dr Anr^ P{r) 



n — 5 



3 (GM^ , MP A , ^ , 
(n + 1) -\^{n- 2)A7irlP, 



(30) 



where the subscript (e) means a quantity evaluated at the boundary rg. Then, 
the dimensionless quantity A is introduced and expressed in terms of homology 
invariants at the boundary r = r^: 



X = 



Er, 



1 



GM"^ 



n — 5 



l{i-("^^^ 



,„-2)| 



(31) 



Figure 2 shows the quantity A as a function of the ratio of the central density 
to that at the boundary, pd Pe- In each panel, lines indicate a series of local 
extrema with a different polytrope index n. Note that all solutions satisfying 
bSq = must lie on this curve. 

At a first glance of Figure 2, we notice that A-curves are bounded from above 
and have peaks in the case of n > 5 {right panel). We call these peaks critical 
points. On the other hand, curves for n < 5 monotonically increase {left panel). 
It follows that, in the case of n > 5, several extremum states correspond to 
a single value of A in some range of A. Thus we can deduce the existence of 
unstable state for n > 5 as follows. Suppose that a polytropic gaseous sphere 
with small radius r^ -^ 0, {pd Pe -^ 1) is stable, i.e., 5'^Sq < (see Appendix C 
for rigorous proof). Due to the continuity, entropy Sq must be local maximum 



10 



along each A-curve. This is true as long as A increases monotonically. When 
the ratio pd Pe exceeds the critical point, however, the extremum configuration 
cannot be stable. If this is stable, the entropy has to be a local maximum 
there, which implies that there exists a local minimum between these two local 
maxima. This contradicts with the fact that all extremum state must lie on 
the A-curve. Therefore, local maximum of the entropy can only be attained at 
the density pd Pe below the critical point and the stability/instability criterion 
is thus obtained from the condition that the quantity 

d\ , , 

= (32) 



ApcIp^ 



first vanishes as increasing the ratio pd Pe- From the monotonicity of p(r), this 
condition implies d\/d^e = 0. Therefore, using the dimensionless quantities 
(22) and the equation (23), the stability/instability criterion reduces to 

= -^ = ^~^ 9i'^e,Ve) .g2^ 

d^e n-5 2Ve ^e 



where 



g[Ue,Ve) 

= 4.f + 2„.., - {8 + 3 (;i±i) } .. - ^ «. + 3 (;i±i) . (34) 

As will be shown in the next section, especially in eq.(53), this criterion can 
be derived also by means of an explicit method, i.e., evaluation of the second 
order variation of the entropy. 

In Table 1, numerical values of dimensionless quantities A and Pc/Pe evaluated 
at the critical point are summarized. Table 1 shows that both two quantities 
decrease as n increases. In the limit n -^ oo(or q ^ 1), they asymptotically 
approach the well-known results of the isothermal sphere [1-3] [5]. 



4 Stability criterion from the second variation of entropy 



In this section, we re-consider the stability/instability of static self-gravitating 
system by evaluating the second variation of entropy. Then, the criterion (33) 
is re-derived, independently of the turning-point analysis. 

According to the principle of maximum entropy, the equilibrium state can be 
attained only when the second variation of entropy 6'^Sq around the extremum 



11 



solution becomes negative, S'^Sq < 0. Conversely, the solution becomes unsta- 
ble if one obtains S'^Sg > 0. The condition 6'^Sq = corresponds to the neutral 
case, in which the equilibrium configuration becomes neither stable nor unsta- 
ble. Thus, the stability/instability criterion can be extracted from S'^Sg = 0. 

Consider the variation of entropy around the extremum configuration, fixing 
the energy E and the mass M. We deal with the density perturbation 6p 
around the polytropic gaseous state p{r). To be specific, we focus on the 
radial mode of the density perturbations under the condition: 



'e 

/ dr ATTr'^5p{r) = 0, 



(35) 



so as to satisfy the mass conservation, SM = 0. In appendix B, variation of 
entropy around the equilibrium configuration is computed up to the second 
order. The resultant expression for 6'^Sq then becomes 



S'S. 



Br, 



+ 



n + 1 

~rp{3/2)/{n-3/2) ^_3/2 

n{n + 1) 1 



^3^pp^$ n-3/2 {5pf 



(„-W3rwi7'''K*"5^^n^ 



(36) 



where the variable (5$ denotes the potential related to the perturbation 5p: 



1 d 

^2 ^^ 



,d6^' 
dr 



ATxG6p{r) 



(37) 



and the quantities -B„ and W are respectively given by 



5. 



5(3/2, n+ 1/2) 



{5(3/2, n-l/2)}("~'/')/("-='/') 
W= fd'^xp'^'/'', 



4V2n 



3.3/2- 

n 

27 



-l/(n-3/2) 



which reduce to i?„ -^ 1 and W -^ M, in the limit of n ^ oo(or q -^ 1). 

Based on the result (36), we now discuss the existence or the absence of per- 
turbation mode satisfying 6'^Sq = for given background p{r). To simplify the 
analysis, it is convenient to introduce the new variable Q{r) defined by 



Sp{r) 



1 dQjr) 
Airr'^ dr 



(3J 



12 



Then the mass conservation (35) imphes the following boundary condition: 
Q(0) = Q{r,) = 0. (39) 

In terms of Q{r), equation (36) becomes 



S'S. 



B„ 



n+ 1 



2^(3/2)/(n-3/2) ^ _ 3/2 



A[Ql 



(40) 



where 



A[Q] 



''^ {5^dQ_ n-3/2 1 (dgV^ 

2Trf7^ n 87rr2pi-i/« \dr ) 



1 n(n + 1) 
?,T^W {n - 3/2)2 



*|-l(^)-^i§ 



■ (41) 



Since both Bn and T are positive, the condition 5'^Sq = is translated to 
A[Q\ = 0. Thus, apart from the pre-factor, we pay attention to the term A[Q]. 
Integrating (41) by parts and using (37), we obtain 



A[Q] --JdrQir) |— + ^^^^^ (4^^2^1-1/n^^ 



Q{r) 



1 n{n + 1) 
3T^W (n - 3/2)2 



rf$ , 3^n-3/2\^^_i_,i/„rfp|^ 



rfr <! -— + 
1 dr 2 



n^ 



(ir 



rfny dr2 Q{ri)K{ri,r2)Q{r2), 




(42) 



where the kernel K is given by 



^■'-^'-.-'--'{^^(iWM^i 



+ 



1 n{n + 1) 
3T2W^ (n - 3/2)2 



F(ri)F(r2) 



with quantity F{r) being 



-I 

(43) 






(44) 
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Therefore, stability of the background configuration p{r) is deduced from the 
following eigenvalue equation: 



dr' k{r,r')Q{r') = KQ{r). 



(45) 



From (45), the stable condition 6'^Sq < indicates that minimum eigenvalue 
Kmin should be positive, while the local minimum of the entropy 6'^Sq > 
means /tmin < 0. Hence, the boundary between stable and the unstable con- 
figuration corresponds to the condition /tmin = 0. Notice that for a given 
polytrope index n, the minimum eigenvalue Kmin of the system with {E, M) 
depends on the radius of the rigid sphere r^. For sufficiently small re, it can 
be shown that the configuration should be stable and we obtain the positive 
eigenvalue, /tmin > (see appendix C). Therefore, to seek the boundary of the 
stability, it is sufficient to investigate the condition /tmin = when increas- 
ing Te. That is, the stability/instability criterion can be extracted from the 
following zero-eigenvalue equation: 



dr' k{r;r')Q{r') = 0, 



which yields 



LQ{r)^ 



d { 1 (P\ d\ n G 

dr \ Airr'^p \ p J dr j n+lr^ 

2 n - 3/2 1 Gm{r 

" 3 n + 1 ~ 

dr'Anr'^Pir') 



Q{r) 



dr'^^^Qir'). (46) 



r'2 



In deriving the above expression, we have utilized the relation (17) and the 
hydrostatic equilibrium conditions, (20) and (21). 

The zero-eigenvalue equation (46) is integro-differential equation which seems 
intractable at a first glance, however, the solution satisfying boundary condi- 
tion (39) is luckily obtained from the knowledge of background configuration, 
p{r) or $(r). To construct the solution Q{r), first note the action of the oper- 
ator L on Anr^'p and m{r): 



L (47rr^p) 



dr 



n Gm{r) P 



n+ 1 



n 



n , ^ n — 3 Gmir) 

-AnGrp ^ ' 



n + 1 
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and 



, , d ( P\ n Gm(r) n — 1 Gm(r) 

Lm{r) = — - + —— ^ = —— ^, 

ar \ p J n + 1 r'^ n+1 r^ 

where we have used the following relations: 

P dlnp n Gm d ( P\ 1 Gm 

p dhir n + lr' dr \p I n + lr^ 



(47) 



From the indication of these equations, we put the ansatz of the solution, 
Q{r) = ci 47rr^p(r) + C2 m^r), (4J 



and determine the coefficients Ci and C2 by substituting (48) into (46). We 
then have 



Te 



, Gm{r' 



Idr'^Q(r') 



n — 3 n — 1 2 n — 3/2 i * /.^x 

^rr ^1 + ^rr ^2 = ^ ^r; = a. 49 

n + 1 n + 1 3n + l / 

rfrMTrr'^Pfr') 







Further, recall that the ansatz (48) must satisfy the boundary condition (39). 
Since the condition Q{0) = is automatically satisfied, the remaining condi- 
tion Q{re) = requires 



47r rfpe ci+ Mc2 = 0. (50) 



Equations (49) and (50) specify the coefficients, which can be expressed in 
terms of the homology invariants (see definitions [25] [26]): 

(n + l)A (n + 1) Me A , ^, 

^1 = 5 '( T^' ^2 = -, -^. (51) 

n — 6 — [n — l)Ue n — 6 — [n — l)Ue 



The solution (48) with (51) seems still uncertain because of the quantity A in 
coefficients, which implicitly depends on the solution (48) itself. This means 
that the value of the coefficients Ci and C2 should be further constrained so 
that the local expression (48) indeed satisfies the non-local equations. In other 
words, the consistency between (48) and A puts the condition for the back- 
ground quantities, pe, Pe and M evaluated at r = r^. This is just the stabil- 
ity/instability criterion we wish to clarify. 
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Now, substituting the solution (48), we evaluate the quantity A explicitly. 



2 n - 3/2 
A = — 

3 n + 1 



Ci / dr — — Aixr'^p + C2 / 



dr 



Gm^ 



(52) 



dr Airr P 



The integrals in the right hand side of (52) can be rewritten repeating the inte- 
gration by parts as shown in Appendix A. Together with the coefficients (51), 
substitution of (55)-(57) into (52) leads to the stability/instability criterion. 
In terms of u-v variables, this gives 



i = H 



n-3/2 



3 n — 3 — (n — l)ue 



{n - 5)ue 



2Ue + Ve 



+ 3 



-Up 



- n — 1 

(n-5){2Me + ^^e} 



2Ue + Ve 



n 



1 



After some algebra, the above equation finally reduces to the following quadratic 
form: 



= 4Mg + 2UeVe 
= 9{Ue,Ve), 



8 + 3 



n 



n 



Up 



n 



3 ^nfn + 1 

-Ve + 3 



n 



(53) 



which exactly matches the criterion derived from the turning-point analysis 

(seeeqs.[33][34]). 

Equation (53) is the main result of our analysis. It determines the critical point 
in the (m, f )-plane, where the local extremum entropy is neither maximum nor 
minimum [5]. 

To see the geometrical meaning of this criterion explicitly, we translate the 
result (53) into the constraint in the {u, f )-plane. In the previous section, we 
see that a family of equilibrium sequences is characterized by the parameters, 
Tg, E and M. These parameters specify the value of A defined in (31). In other 
words, such a polytrope system must lie on the straight line: 



n 



(n-5)A-3/2 



n 



2 3 

u 

n-1 2. 



(54) 



from equation (31). Note also that the polytrope solution must lie on the u- 
V trajectory as shown in Figure 1. Hence, the equilibrium state with fixed 
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parameter set (re, E, M) can exist only if the u-v curve intersects the straight 
hne (54). From the configuration of the u-v trajectory, we notice that the 
quantity A is bounded from above, A < Ac in the case of n > 5. That is, for 
the fines (54) witfi A larger tfian Ac, tfie intersection ceases to exist. Since 
tfie critical value Ac must satisfy d\/dr = 0, tfiis yields (33) or (53) and tfie 
condition g{ue,Ve) = determines tfie critical point along eacfi u-v trajectory. 

In Figure 3, witfi a great advantage of m, f-variables, tfie stability/instability 
criterion (53) is examined. In eacfi panel, tfie tfiick-solid lines denote tfie tra- 
jectories of tfie Emden solution. Tfie tfiin-solid lines correspond to tfie equation 
g{u,v) = 0. From Figure 3, one can observe tfiat tfie curve g{u,v) = dras- 
tically cfianges its befiavior around tfie index n = 5 (compare n = 4.8 witfi 
n = 5.2 case). Note also tfie different asymptotic befiavior of tfie Emden so- 
lutions between tfie n > 5 cases and tfie n < 5 cases (see also Figure 1). As a 
consequence, tfie critical point appears wfien tfie polytrope index n > 5 and 
it disappears in tfie n < 5 cases. 

For tfie n > 5 cases, we also plot tfie straigfit line (54) witfi tfie critical value Ac 
(dasfied line). It sfiould be empfiasized tfiat tfiree lines, (i.e. tfie trajectory of 
tfie Emden solution, tfie stability/instability criterion (53) and tfie condition 
(54) witfi tfie critical value Ac) intersect at tfie critical point. For comparison, 
tfie result of tfie isotfiermal gas spfiere is plotted in Figure 3. Tfie existence of 
tfie critical point still fiolds in tfie limit n — > cx) [1-3] [5]. 



5 Discussion and conclusions 



In tfiis paper, we fiave applied Tsallis-type generalized entropy Sg to tfie prob- 
lem of seeking tfie stable distribution of self-gravitating systems. In contrast 
to previous work using tfie Boltzmann-Gibbs entropy, tfie local extremum 
state of Tsallis-type entropy fias been found to be equivalent to tfie stellar 
polytrope system and Tsallis' g-parameter is related to tfie polytrope index 
n (see eq.[14]). Tfien we apply tfie usual turning-point analysis to explore 
stability/instability of tfie local extremum state. We note tfiat tfie fiomology 
invariants {u, v) are very useful to investigate tfie nature of equilibrium state 
of tfie self-gravitating system, including tfie stability/instability criterion (33). 
Tfien we develop tfie second variation of entropy and re-examine tfie stabil- 
ity/instability of tfie local extremum state. Tfie stability/instability criterion 
(53), wfiicfi exactly coincides witfi result (33), can be obtained by solving tfie 
integro-differential equation of zero-eigenvalue state. 

Tfie results imply tfie important conclusion tfiat tfie polytropic gaseous spfiere 
witfiin radius Vg exfiibits tfie gravotfiermal instability in tfie case of polytrope 
index n > 5. Tfie cfiaracteristic values A and Pc/Pe are also evaluated at 



17 



critical point and presented in Table 1. These values asymptotically approach 
the well-known result of the isothermal sphere in the limit n -^ oo(or q — > 
1). In the isothermal case, the gravothermal instability can be interpreted in 
terms of a negative specific heat which seems to be a common feature in the 
self-gravitating systems[2]. We note, however, that the system should have 
sufficient amount of outer normal part in order to trigger the instability; the 
specific heat for total system should be positive for onset of the instability, 
although its central part has negative specific heat. A similar argument holds 
for the polytropic system so that the instability appears for the system of 
n > 5 which can have an elongated outer part. Evaluation of the specific heat 
based on the canonical ensemble corresponding to Tsallis' type entropy seems 
to be very important on this issue[31]. 

Finally, we comment on the role of constraints M and E to the extremum solu- 
tion (9). In our present analysis, standard definition of mean values is adopted 
for mass and energy (see eqs.[3][4]). Tsallis, Mendes and Plastino[17] recently 
showed that this choice yields undesirable divergence in some physical sys- 
tems including Levy random walk. To overcome the mathematical difficulty, 
they suggest that the g-generalized mean value should be used with a correct 
normalization instead of the standard mean value (see also [18]). According 
to their suggestion, in Ref . [28] , g-generalized, but un-normalized mean value 
was used to define E and M . In this case, the resultant form of extremum 
state essentially remains unchanged and the polytropic equation of state still 
holds. We note, however, that if we use the normalized g-expectation value 
correctly, the problematic difficulty conversely arises. This is due to the fact 
that the potential energy U associated with the distribution function /(a;, v) 
(see eqs.[4][5]) becomes nonlinear function of /, while the new formalism in 
Refs.[17,18] implicitly assumes the linear function. Certainly, the normalized 
g-expectation value will play an essential role to avoid the unexpected singu- 
lar behavior, however, no undesirable divergence has appeared in our present 
analysis. Furthermore, it has been shown that the Tsallis formalism with stan- 
dard linear mean values still verifies the Legendre transform structure, leading 
to the standard results of thermodynamic relation [15,16]. Therefore, at least 
in our case of the self-gravitating system described by the distribution func- 
tion /, the analysis using standard mean values can be validated and the 
conclusion remains correct. 

The framework of nonextensive statistics based on Tsallis' type entropy seems 
to give a consistent generalization of the usual thermodynamical structure 
[14,17-19,23]. These works, however, are mainly concerned with construction 
of a consistent formal framework. In the light of this, our specific application 
of the Tsallis-type entropy to the self-gravitating systems is the first realistic 
physical consideration. Due to the spherical symmetry, the system is sim- 
ple and we can easily evaluate physical quantities, e.g., pressure, energy and 
entropy. Nevertheless, the result still shows a very interesting phenomenon. 



i.e., the gravothermal instability. Furthermore a fact of the gravity being a 
long-range force suggests that the self-gravitating system is one of the most 
preferable and interesting testing grounds for the framework of nonextensive 
statistics. 
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Appendix A: Some formulae for integration 

Here we list some formulae which have been used in section 3 and 4. From the 
hydrostatic equations (20), (21) and the integration by parts, we obtain 



/Gtti f dP f 

dr -^ A-nr^p = - dr Airr^ — = -AnrlPe + 3 dr Anr"^ P, (55) 



Gm^ GM^ „?, . ^Gm GM^ „?, , ^dP 



/ dr — r— = \-2 / drAnr^ p = 2 / drAnr'^-— 

J r^ Te J r Te J dr 



- Svrr^Pg + 6 f dr Anr^ P, (56) 







Similarly we can rewrite the integration of pressure P as 



, A 2 n 47r o ^ f ^ An ^dP 47r o ^ G ( dr d , o. 

drAirr^ P= — r^l P^ - dr — r'^ — - = — r-^ p -\ / -(m^^J 

3 "" J 3 dr 3 ^ Q J r dr^ ' 



47r o „ GM2 n + 1 ? , d ( P\ 

= Y'^'^^^^^-^j/'d^\-^)'^ 

An 3 ^ GM^ n + 1 [ MP, . ^ . 
-rt P, H < - / dr Anr^ P 



3 ' ^ Qr, 6 1 p. ^ 



which leads to 



drAnr^P = — UnrlP, - (n+ 1)^^ + ^| . (57) 

n - 5 I pe re I 



Appendix B: Second variation of entropy 



In this appendix, we derive the second variation of entropy around the equi- 
librium state 6Sg = 0. 

We first express the entropy (2) of the extremum state. Substitution of the 
distribution function (16) leads to 
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C<(max) 



Q 



1^ I rf^^rf^t; (/«-/) 



3 

n 

2 



M 



Br, 



2^(3/2)/(n-3/2) 



/ Sxp 



1+1/n 



(5J 



where the constant Bn is given by 

5(3/2, n+ 1/2) 



B„ 



{5(3/2, n- 1/2)} 



(n-l/2)/{n-3/2) 



. o. 3/21-l/{n-3/2) 

4V27r U - - j 



.(59) 



Just for convenience, we introduce the quantity W 



(60) 



and vary (58) with respect to p and T keeping E and M fixed. Up to the 
second order, we obtain 



^^-^ 



Br, 



T 



_3 ri^(5T + 5T51^ 

2^(3/2)/(n-3/2) ^2 j 



+ [n--] 5W 



n 



W 



^n-3/2; 2T2 



(<5T)^ 



(61) 



On the other hand, the constraint 5E = gives 



"-^^-^K^l^'^'^^/^^-"* 



3 /n-3/2 

2 I n + 



^j 5{TW)+ f £x(^6p+-SpS^Y 



(62) 



We thus find 



ST (W + SW) = - 



2 / n+ 1 



3 In- 3/2^ 



/ d^x Usp + ^SpS<^) - TSW. 



(63) 



(In arriving at eq.[62], we have used the fact that J d'^x ^6p = J d'^x p6^). 
Substituting (63) into (61), we ehminate the quantity 6T. Then, collecting the 
second order terms yields the second variation of entropy S'^Sq as follows: 



s's., 



B„ 



2^(3/2)/(n-3/2) 



n + 1 

n-3/2 



,3 6p 5$ 
2T 



+ n6W 
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3 1 n (2 n+1 



/ d^x <^5p + T5W 



(64) 



4 T^W n- 3/2 [3 n - 3/2 
Now consider the variation of W: 

5W = 6^j £x pi+i/") = ^ / d'x pV« Up + ^ ^) , (65) 



Neglecting the higher order contributions and keeping the second order terms 
only, substitution of (65) finally leads to the expression (36): 



s^s, ""- " + ' 



9 2^{3/2)/{n-3/2) n-3/2 

n{n + 1) 



J5p5^ n-3/2_i6pl 
\ 2T ^ 2n pi-iA 



+ 



;n- 3/2)2 32-2p^ 



L,{/..(.,|^.,.)^- 



. (66) 



Note that the above equation indeed reduces to the well-known result in the 
isothermal sphere. Using the fact that B^ -^ 1 and W -^ M in the limit 
n — * oo, we obtain 



525/-^ - £x{ -^ + ^^^^ - — ^ / £x^5p 



which is equivalent to the expression (16) of Antonov[l] and equation (ALIO) 
of Padmanabhan[5]. 



Appendix C: Positivity of minimum eigenvalue k^. 



For the sake of completeness, in this appendix, we will prove that for suffi- 
ciently small Tg, minimum eigenvalue of the equation (45), /tmin can become 
positive, irrespective of the polytrope index n. 

First recall that from (42) and (45), positivity K,„iin > is equivalent to the 
condition A[Q] < 0. This gives 



''' , l6<^dQ n-3/2 1 ( dQ^ ^" 

dr { — z-— — h 



2T dr n Sirr'^p^^^/'^^ \ dr 

' r \ 
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n[n + 1) 



'iT'^W (n - 3/2)=^ 




3|^n-3/2\ ^,/ 1 rfQ 
n I \ dr 



< 0. (67) 



In the above expression, while the second term in the right hand side of equa- 
tion is always negative, the first term can be expressed as 

[1st term] = -±^(H-l) J dr ^, (68) 



where we define 




(69) 



Thus, the inequality H > 1 provides a sufficient condition for the positivity 

OI /tmiii' 

Now, we rewrite the definition (69) integrating by parts. We obtain 

_^( 1 (£)^Uff^^. (TO) 

dr I Anr'^p \p I dr \ n+1 r"^ 



The above equation can be regarded as the eigenvalue equation with eigenvalue 
H. In fact, one can show that equation (70) has minimum eigenvalue -ffmin = 1 
for some small radius rg. Consider the function 

Qmin = c |47rrV(r) —m{r) \ = c (u — - j m, (71) 



where c is merely a constant. The above function behaves like Qmin — c (2n/(n— 
l))m ^ near the origin r = and it also becomes vanishing at r = ri, where 
u{ri) = {n — 3)/(n — 1). Further, Qmin does not vanish during the interval 
(0,ri). Hence, when setting the radius of the wall as rg = ri, the function 
(71) indeed corresponds to the ground state of the eigenvalue equation, i.e, 
perturbation mode without any nodes. Therefore, this fact proves that if we 
suitably choose a smaller radius r^ < ri, the eigenvalue of the system (70) 
should be larger than unity, i.e, H > 1. Hence, all equilibrium configurations 
with Tg < ri become stable. 
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Table 1 

Energy-radius-mass relation and density contrast between center and edge evaluated 

at a critical point for given polytrope index n or q 



n 


q 


/ 
A 
V 




N 
r 
/ 


Pe 


5 


9 

7 










6 


1.22 




3.78 




1.27x10^ 


7 


1.18 




1.69 




1.46x10^ 


8 


1.15 




1.162 




4.83x10'' 


9 


1.13 




0.932 




2.40x10^ 


10 


1.12 




0.804 




1.46x10'' 


30 


1.04 




0.429 




1590 


50 


1.02 




0.388 




1130 


100 


1.01 




0.360 




887 


CX) 


1 




0.335 




709 
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Fig. 1. Trajectories of Emden solution in {u,v)-p\ane 
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Fig. 2. Energy-radius-density contrast relationship 
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Fig. 3. Stability/instability criterion in (u, t;)-plane 
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